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The possibility of the existence of states with a spontaneous interlayer phase coherence in multi- 
layer electron systems in a high perpendicular to the layers magnetic field is investigated. It is shown 
that phase coherence can be established in such systems only within individual pairs of adjacent 
layers, while such coherence does not exist between layers of different pairs. The conditions for 
stability of the state with interlayer phase coherence against transition to a charge-ordered state are 
determined. It is shown that in the system with the number of layers TV < 10 these conditions are 
satisfied at any value of the interlayer distance d. For N > 10 there are two intervals of stability: 
at sufficiently large and at sufficiently small d. For N — > oo the stability interval in the region of 
small d vanishes. 

PACS numbers: 73.45.Cd, 73.21.Ac, 73.21. Cd 

I. INTRODUCTION 

In the last 10 years main attention in studying the quantum Hall effect and related phenomena was given to bilayer 
systems. One of the most interesting features of bilayer quantum Hall systems is the possibility of realizing in them a 
superfluid state of a gas of electron-hole pairs with spatially separated components. For the first time, the possibility 
of superfluidity of this type was predicted in Refs. y] and with ref erence to bilayer electron-hole systems in zero 
magnetic field (see also Refs. and0). In Ref. E (see also Refs. IfiFHril for more details) it was shown that strong 
magnetic field applied to a bilayer electron-hole system favors the formation of a superfluid state, since, unlike the 
situation considered in Refs. QandH in the latter case it is not necessary to satisfy the condition of nesting of the 
Fermi surfaces of the electrons and holes. 

Owing to the electron-hole symmetry, a bilayer electron-electron system in a quantizing magnetic field with filling 
factors of the layers equal to v\ = v and = 1 — v [y < 1) is a practically complete analog of bilayer electron-hole 
systems, and a transition to a superfluid phase is possible in such a system as we lii2iiiii2ii&i2ii£ii&. Recently, interest 
in this question has risen significantly in connection with experimentsi^i^ in which an effect of the Josephson type 
was observed (a sharp peak of the differential tunneling conductance at gate voltage biased near zero). Furthermore, 
in a recent experiment — on the interlayer drag in bilayer systems a strong suppression of the longitudinal drag current 
and an appearance of the Hall component of the drag current was observed. This effect takes place in the situation 
when the total filling factor v to t — v\ + v 2 becomes close to unity. The effects observed in Refs. d 18lll9l can be 
explained on the assumption that a spontaneous interlayer phase coherence arises in the system and a transition of the 
system to a superfluid state occurs. A direct experimental proof of Bose-Einstein condensation of interwell excitons 
(photoexcited electron-hole pairs with spatially separated components) in double quantum wells in zero magnetic field 
was also obtained recently^. The prooSS was based on the measured temperature dependence of the luminescence 
• i-h , spectra. 

For further understanding of this phenomenon it is naturally to study the superfluid properties of multilayer electron 
systems. Multilayer systems are of particular interest for studying collective properties of electron-hole superfluids 
by optical and acoustic methods, since the integrated intensity of the interaction of such systems with external fields 
increases in proportion to the number of layers. The present study is devoted to investigation of the effect of interlayer 
phase coherence in quantum Hall superlattices. 

The realization of multilayer quantum Hall systems was reported in Refs. The research3ii22i2£ was mainly 

concerned with analysis of the properties of chiral edge states in such systems. Previously, the question of interlayer 
phase coherence in multilayer quantum Hall systems was considered in Refs. 00] and|25[ The authors of Ref. [25| have 
shown that a state in which phase coherence arises between all layers can be realized in such systems (below, this is 
referred to as a state with global phase coherence). It has been shown in Ref. that in a multilayer quantum Hall 
system another state can arise (which we shall call a dimer state) in which the system separates into pairs of adjacent 
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layers and the interlayer phase coherence is established only within each pair, while there is no such coherence between 
layers belonging to different pairs. This raises the question of reconciliation of the results of Refs. and|2^. In the 
present paper we show that the states of Refs. and|2^ correspond to two different solutions of the self-consistency 
equations for the order parameters describing the interlayer phase coherence. We find that the energy of the dimer 
state is less than the energy of the state found in Ref. |25| and, hence, it is the dimer state that should be considered 
as a candidate for the grou nd state of such a system. 

It was shown in Ref. |25| that conditions for the existence of the state with interlayer phase coherence in multilayer 
systems are more restrictive than for bilayer systems. Multilayer quantum Hall systems, unlike the bilayer ones, 
demonstrate a tendency toward a formation of a charge-ordered state (where the filling factors of adjacent layers 
become equal to and 1). In Ref. [25|the stability conditions for the state with interlayer phase coherence were analyzed 
under the assumption that the phase coherence has a global character. Since this assumption is not confirmed, the 
stability conditions will have to be reexamined. We find that in the limit of an infinite number of layers (the case 
considered in Ref. I25I) the dimer state is stable against transition to the charge-ordered state already at interlayer 
distances d > d c as 1A5£b (£b is the magnetic length), which is less than the critical value (d c ~ 1-7£b ) obtained in 
Ref. E 

In connection with the fact that the properties of systems with number of layers N — 2 and N — > 00 are substantially 
different, and in real physical systems the number of layers is always finite, an important question arises as to how 
the stability conditions depend on N. In this paper this question is studied with reference to the system with even N. 
It is found that for N < 10 the state with interlayer phase coherence is stable for arbitrary d. For N > 10 stability 
takes place only for d < d c \ or for d > d C 2, where d c \ decreases and d C 2 increases with increasing N while remaining 
in the interval < d c i < d c2 < d c . 



II. INTERLAYER PHASE COHERENCE IN BILAYER SYSTEM WITH TOTAL FILLING FACTOR 

Vtot = 1 

Prior studying the multilayer system, we describe the approach used in this paper with reference to the simplest 
case N = 2 and recall the bilayer situation. Let us consider a double-layer electron system subjected by a high 
magnetic field applied perpendicular to the layers. We imply the filling factor v to t — v\ + ^2 = 1. Tunneling between 
layers is assumed to be so weak to be neglected in the Hamiltonian. If the cyclotron frequency is much higher than all 
the other characteristic energies of the problem, the system can be treated in the lowest Landau level approximation, 
in which the influence of the upper unfilled Landau levels on the dynamics of the system is neglected. 

The Hamiltonian of this system has the form 

1 
2 

n,n' = l,2 



H = \ I /rf 2 ^V„^+(r n )*+(r„OKn'(|rn-r„H)*„'(rn')*n(rr l ) 

- / d 2 r n ^+{r n )m n (v n ) +H BG , (1) 

n=l,2 ^ 

where ^ n are the Fermi field operators, 

e 2 

V m '(r) = m - (2) 

e [r 2 + (n — n') 2 d 2 \ 

is the Coulomb potential, e, the dielectric constant, and /i„, the chemical potential in layer n. The term Hbg in Q 
describes the interaction of electrons with positively charged impurities in doping layers. For simplicity we assume 
that the z coordinates of the electron layers and doping layers coincide. If one takes into account the difference of 
the z coordinates of the electron and doping layers which takes place in the real physical situation it results in an 
appearance of additional constant term in the energy that does not influence the effects studied in this paper. 

For further analysis it is convenient to transform to the operators of the Fourier components of the electron density 
/3„(q), in terms of which the Hamiltonian Q can be written in the form 



n.n' = l,2 q 



q 2 £ 2 

Pn(qW(-q) - <W exp( 2^)/?n(0) 



VnPn(0)+H B G, (3) 
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where S is the area of the layer. The operator p n (q) expressed in terms of the creation and annihilation operators for 
electrons in the lowest Landau level (a + (X), a(X)) (X is the guiding center of the orbit) has the form: 



,(q) = J2 <( X + g«4/2)a»(* - 9/1/2) exp(iq x X - q 2 £ 2 B /4), 



x 



In Eq. © the quantity 



V nn > {q) 



2ire 2 
eq 



cxp(— qd\n — n'\) 



(4) 



(5) 



is the Fourier component of the Coulomb potential. 

For further analysis we use the mean field approximation. In this approximation the Hamiltonian @ takes the 
form 



Hmf = {^iat(X)ai{X) + e 2 a+(X)a 2 {X) - [JAa+ (X)a 2 {X) + h.c.}} 
The energies e n in © arc 



x 



Vn - IVr, 



Mr, 



(6) 



(7) 



where V n are the Fourier components of the Coulomb interaction at q = screened by positively charged doping 
impurities. For the bilayer system these quantities are equal to 



V x = 
V 2 = 



^ — -pr h m 

2tt£ 2 b q^O 



In (7) [vx-\ 



V 22 {q) v 2 



1 i 2 ( q ! ( V 2 - - 



V21 (q) v\ 



= Wv, 



(8) 



The parameter v = (y\ — v 2 )/2 describes the value of the imbalance of the filling factors of the layers. In Eq. © A 
is the order parameter, defined as A = (a 2 (X)ai(X)} . In general case the order parameter A = | A|e ztp is complex 
valued. We restrict consideration to the case where its modulus |A| and phase ip arc independent of X. In deriving 
Eq. JBJ) we have also taken into account that the averages (a„(X)a n (X)) are equal to v n . In Eqs. iJBJ, Q, and (JSJ 
we have used the following energy parameters: the parameter 



W 



that describes the energy of the direct Coulomb interaction between layers, and the parameters 



/=i S y 1 , M exp(V4/2) = \/f£, 



(9) 



(10) 



J - |E^( 9 )exp(-^/2) = ^exp (£) e rfc (^) (11, 

that describe the energy of the intralayer and interlayer exchange interactions. 

In a state in which the order parameter A is nonzero, the motion of electrons in one layer is correlated with the 
motion of holes (empty single-particle states on the lowest Landau level) in the other layer. This state can be treated 
as a gas of electron-hole pairs, which are composite bosons. The Bose-Einstein condensate of such pairs (a true 
condensate at T = or a quasi-condensate with a fluctuating phase at T ^ 0) will demonstrate superfluid properties. 
The state with nonzero superflow is the state in which nondissipative electrical current in one layer is accompanied 
by equal and oppositely directed nondissipative electrical current in the other layer. 

For finding the self-consistency conditions for the order parameter A and chemical potentials fix and /i 2 we subject 
the operators a\ and a 2 to a u — v transformation of the form 



a\ = ua + v*(3 + , 
a 2 = u*/3 + — va. 



(12) 
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The operators a and j3 will satisfy the Fermi commutation relations if |u| 2 + \v\ 2 = 1. This allows to seek the 
coefficients of the u — v transformation in the form u — cos(0/2) and v — sin(0/2)e lx . From the requirement of 
vanishing the terms nondiagonal in a and f3 in the transformed Hamiltonian, we find that 

s inP)= — ^ |A| rn*P> = — 6 (13) 

and that x coincides with the phase tp of the order parameter. Here we introduce the notation e = (ei — o^)/2. 
Using Eqs. I|12fl and we obtain the following self-consistency equations: 

H + v 2 = 1 = (aj-oi) + (a+a 2 ) = 1 + iV F (£ Q ) - N F (£p), (14) 

2£ = (a+ ai ) - (a+a 2 ) = —_=!_— [1 - AT F (£ a ) - A^)], (15) 

V e + ( J I A I) 

A = (atai) = — ; JA [1 - N F {£ a ) - N F (£a)}, (16) 



where 



8 am = #HW ± - (17) 



are the energies of elementary excitations and Np(£) — [exp(£/T) + 1] _1 is a the Fermi distribution function. It 
follows from Eq. i|14|) that £ a — £p , which leads to the condition ei + £2 = 0. From Eqs. (|15fl and Ijltjfl one can easily 
find the dependence of the modulus of the order parameter on the value of the imbalance and temperature. Below 
we specify the case of zero temperatures, when this dependence has the form 



|A(iO| = ^-i>2. (18) 

The differential of the free energy of the system at V\ + z/ 2 = 1 and T = is equal to 

dF = nidvi + ii^dvi = 2jldv^ (19) 

where fl = (/xj — /ia)/2. Here and below we give the energy per electron. If the analytical form of the function jl{v) 
is known, the energy of the system as a function of the imbalance value is easily found: 

F(v) = F(-l/2) + 2 / di>'p,(i>'), (20) 
J -1/2 

where F(— 1/2) is the energy of the system for v = —1/2, i.e., for v\ = and = 1- At such filling factors the energy 
of the system is the sum of the energies of the direct Coulomb interaction between layers and the intralayer exchange 
interaction in layer 2, i.e., F(—l/2) — W/A — 1/2. Using Eq. Q and taking into account the relation e = — Jv that 
follows from Eqs. (|15f) and (|lt>|) . we find the following expression for fr. 

fi(p) = {W-I + J)v. 

Accordingly, 

F(v) = - I -±^ + (W-I + J)v 2 . (21) 
It is useful to present expression (|21|l in the form 

F = WD 2 -^[vl + il-v,) 2 ]- J|A| 2 (22) 

from which the physical meaning of each of the terms becomes obvious. The first term is the energy of the direct 
Coulomb interaction at filling factor imbalance v. The second term is the sum of the energies of the intralayer exchange 
interaction in layers 1 and 2. The third term is the energy connected with the interlayer phase coherence. 
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Since energy (1221 is independent of the phase ip of the order parameter, the interlayer phase coherence is spontaneous 
one. The minimum energy of a state with interlayer phase coherence is reached for the modulus and phase of the 
order parameter independent of the coordinate. A state with nonzero superfluid current that corresponds to nonzero 
gradients of the phase is higher in energy, but in the case where the gradient of the phase is less than certain critical 
value (which depends on the ratio d/is), such a state remains metastable, because the system cannot emit normal 
excitations, i.e., the Landau criterion of superfluidity is satisfied. 

To conclude this Section we turn our attention to an important observation that follows from formula (|21|l . Using 
Eqs. lj§|l- i|ll|) . we see that the coefficient of the v 2 term on the right-hand side of (|21|> is positive for arbitrary d ^ 0. 
This plays a fundamental role for the existence of the state with interlayer phase coherence in bilayer systems. Indeed, 
if the sign of the expression W — I + J were negative, then the system would go to a state with a maximum imbalance 
of the filling factors, v = ±1/2, for which |A| = (see Eq. (|T%10 and, consequently, interlayer phase coherence is 
absent. Skipping ahead, we note that in multilayer systems, in contrast, such a mechanism of destroying the phase 
coherence can indeed be realized in a certain range of d. This question is analyzed in details in Sec. IIVI 



'm = */~«p(^ (24) 



III. STATES WITH INTERLAYER PHASE COHERENCE IN MULTILAYER SYSTEM 

Let us consider a system with the number of layers N — ► oo with an average filling factor per layer equal 1/2. The 
Hamiltonian of such a system can be written in the form @ with the sum over two layers replaced by a sum over an 
infinite number of layers. In the mean field approximation, the Hamiltonian of the multilayer systems is reduced to 

Hmf = ^ X^ e " a ™ a ™ _ \ X! [^™A(n, m)a+a n+m + h.c.]}, (23) 
where the energy parameters J m are equal to 

erfc (^] 

(the parameter J introduced in the previous Section coincides with Ji) and the order parameters are defined as 

A(n, m) = (<z+ +m a n ) (25) 

In general case the quantities A(n,m) can be nonzero for arbitrary m (not only for m = ±1). This corresponds to 
the situation when a correlation arises between the electrons belonging not only to the adjacent layers but also to 
the layers arbitrarily far apart. Since the Coulomb potential is long-ranged, such correlations may give significant 
contribution to the energy, and it is important to take them into account for finding the ground state of the system. 

It is natural to expect from physical reasons that if an imbalance of the filling factors arise in such a system it has 
a periodic character. Let us consider the situation when such a period does not exceed twice the distance between 
adjacent layers. As it shown below, the doubling of the period may happen even in the absence of imbalance of the 
filling factors. Latter possibility corresponds to the state with interlayer phase coherence of the dimer type. 

Let the filling factors of the layers be equal to v n = 1/2 + (— l) n D. In this case the quantity v can be treated as 
a charge-ordering parameter. Basing on the assumption of periodicity of the system, we seek the order parameter 
A(n,m) in the form 

A(n,m) = Ai(m) + (-l)"A 2 (m). (26) 

In general case the quantities A^^m) can be complex valued, i.e., they can contain phase factors e* Vl < 2 >( m ). The 
phase factors can no longer be specified independently, since the self-consistency equation imposes certain conditions 
on the phase difference <Pi(2)(w). However, these conditions do not fix all of the phases, and a certain arbitrariness 
in the choice of the quantities </?i(2)(m) remains. It can be shown that the energy of the system is independent of the 
phases of the order parameters if the latter satisfy the self-consistency equations (we stress that we are considering 
the case when the modulus and phase of the order parameter are independent of X). An arbitrariness in the choice 
of phases is just a reflection of the spontaneous character of the phase coherence. It allows us to be more specific and 
restrict consideration to the simplest case of all the phases are equal to zero and the quantities Ai(m) and A2(m) are 
real valued. 

Transforming in Eq. (23) to the Fourier components of the operators a n 

a » = 7=E ei92r v> ( 27 ) 
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we obtain the following expression for the Hamiltonian of the system: 

H MF =J2^{[e-MQz)K, a Q* + [Z- A 2r{qz)-i&2i{qz)]ai+^a q J. (28) 

X 

Here the following notations are introduced e = (e2k +£2fc-i)/2, e = (tik — £2fe-i)/2. In view of the assumed periodicity 
in z the quantities e and e are independent of k. The functions Ai(q z ) and A 2r (i)(<7z) in (|28|l are defined as 

oo 

Ai(fc) = 2 ^ J m Ai(m) cos(mg z ), (29) 



A 2r (g z ) = 2 ^ J 2m A 2 (2m) cos(2mg z ), (30) 

m= 1 



A 2l fe) = 2 ^ J 2m+1 A 2 (2m + 1) sin[(2m + l)gj. (31) 



m=0 

Diagonalizing the Hamiltonian (|28|l and calculating the averages of iT-'Ot , we arrive to the self-consistency conditions 
for the order parameters. Here we present these equations for the case T = 0: 

i r' 2 e-AM 

27r J-tt/2 E{q z ) 

Ai(2rn) = 0, (33) 

Al(2m+1) ' r^ ^t^fe) , (34) 

A fOml 1 f /!> W cos(2mg z )[g-A 2r ( gz )] 

A 2 (2m) = - — / --^ , (35) 

A ,„ 1 f /2 , sin[(2m + l) gz ]A 2i ( gz ) 

A 2 2m + 1=— / dq 2 -— . (36) 

In equations (|32|) and (|34[1 - H36|> the spectrum of energies of elementary excitations E(q z ) is determined by the following 
relation: 



E(g z ) = ^[g- A 2r (g 2 )] 2 + A\{q z ) + AUq z ). (37) 

The state with global phase coherence considered in Ref. |2^ corresponds to the solution of i|32[l - H3t)|) with all the 
quantities A2(2m + 1) equal zero (in this case Eq. (|36|l is satisfied automatically) and Ai(2m + 1)^0 for arbitrarily 
large m. The latter means that phase coherence is established between all layers. 

It is easy to see that the self-consistency equations (|32|) -(I36 |I are also satisfied by a solution of an essentially different 
type: Ai(±l) = ±A 2 (±1) ^ 0, but Ai(m) = A 2 (m) = at m ^ ±1. In this case A(2fe, 1) = A(2fc + 1, -1) ^ (or 
A(2fc —1,1) = A(2k, —1) 7^ 0), but all the remaining A(n, m) are zero. In the state corresponding to such a solution, 
the system separates into pairs of layers, and the interlayer phase coherence is established only within each pair. We 
call this state the dimer state. 

Writing the differential of the free energy of a multilayer system and integrating it over i>, we obtain an expression 
for the free energy of the multilayer system: 

F{v) = V c (p) - I (j + i)A - 2 Jm[A\{m) + A 2 2 (m)}. (38) 
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where the first term V c (v) corresponds to the energy of the direct Coulomb interaction, which is independent of the 
concrete form of A^^jm). 

Using the answer 1)38(1 , let us compare the energies of the state with global phase coherence and of the dimer state 
for v — 0. At zero imbalance Eqs. 1(32(1 - 1(361) give the following expressions for the order parameters Ai^m). For a 
state with global coherence we have 

. , , sin(7rm/2) . . . n 

Ai m = — — , A 2 m) = 0. 39 

mir 

For the dimer state Ai(m) = ±A2(m) = 5i ,i m i/4. Consequently, the energy difference of the dimer state and the 
state with global phase coherence satisfies the inequality 

- -t + hp*»wkw ~ "T f 1 - ht <sW) - < 40 > 

k — \ k — / 

In deriving Eq. I|4(J|) we take into account that for d ^ the parameters J m the smaller the large |m|. For d = 0, 
when all J m are identical, the energies of the two states under consideration are equal. A numerical analysis of the 
self-consistency equations shows that the inequality Fj, < F g holds for arbitrary values of v and d ^ 0. Since the 
energy of the state with global phase coherence is higher than the energy of the dimer state, this state but not the 
state considered in Ref. |25j should be considered as a possible candidate for the ground state of the system. 



IV. COMPARISON OF THE ENERGIES OF THE DIMER AND CHARGE-ORDERED STATES 

Let us now consider the question of the stability of the dimer state in a multilayer system with even number of 
layers against a transition to the charge-ordered state. We should specify the explicit form of the energy of the direct 
Coulomb interaction K(i'). It can be easily done by using the fact that this energy is equal to the energy of the 
electric field induced inside the superlattice due to deviation of the filling factors of individual layers from the average 
value. Since the system on the whole remains electrically neutral, the electric field outside the superlattice is equal 
to zero. 

For instance, for the filling factors of the layers equal to z^2fc-i = 1/2 — v and V2k = 1/2 + v (k = 1, 2, . . . , N/2), 
the electric field is nonzero between layers 1 and 2, 3 and 4, 5 and 6, etc., while between layers 2 and 3, 4 and 5, etc. 
it equals to zero. Then, the electric field energy per electron, as in the bilayer system, would be equal to V c = Wv 2 . 
However, in contrast to the bilayer system, in the multilayer system substantial benefit in the energy of the direct 
Coulomb interaction can be achieved by changing the filling factors of only the two outer layers. Indeed, if the filling 
factors of the outer layers take the values v\ = 1/2 — v/2 and vn = 1/2 + P/2, then the electric field is induced between 
all the layers, and its absolute value would be smaller by a factor of two than in the previous case (the direction of 
the field would change sign on crossing a layer). Consequently, in the last case the energy V c (energy per electron) is 
in two times smaller than in the first case (with an accuracy up to 1/N). In such a state the total energy is 



F = - 



I+J 



W 
2 



J\v 2 . 



(41) 



In deriving Eq. (|41|l we take into account that in the dimer state the order parameter A = 2A;l(2) has the same 
dependence on v as in the bilayer system (Eq. 1(18(1 ~). In contrast to Eq. 1(21(1. the sign of the coefficient of the v 2 term 
in Eq. 141(1 depends on d: for d < lAhls this coefficient is negative, while at larger d it is positive. Consequently, for 
d < 1A5£b and N — > oo the minimum energy corresponds to the state with the maximum charge ordering, in which 
all the A(n,m) go to zero and phase coherence is absent. In the opposite case d > 1A5£b the minimum of energy 
((41(1 is reached for v — and charge ordering becomes energetically unfavorable. 

The above analysis is easily generalized to the case of finite even N. Let us consider the configuration where the 
inner layers are in a completely or partially charge-ordered state, and for these layers the imbalance of the filling 
factors is equal to £, n . The imbalance of the filling factors of the outer layers v ex = (vn ~ v i)/2 can be different from 
i>i n . We imply that phase coherence is established only between the inner layers and has the dimer character. The 
energy (per electron) of such a state is equal to 



F{v ex ,v m ) = W 



N 



N 



(pi: 



N 



N 



1 o\ 2 

- + vL 



N V4 



1 



4 



J- 



N -2 (1 



N 



Pi 



(42) 



The minimum of F{y ex , v in ) should be found in the domain specified by inequalities \v in \ < 1/2 and \v ex \ < 1/2. The 
function F(p ex ,i>i n ) has only one extremum (or saddle point) at v ex = £j„ = 0. Therefore, its minimum is reached 
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FIG. 1: Domains of existence of the dimer (D) and charge-ordered (CO) states and d c i ( o), d C 2 (•) as functions of the number 
of layers N. The interlayer distance d is given in units of the magnetic length £b ■ 



either at the point of the extremum or at the edge of the domain of allowed values of £^ n , i> ex . Thus, for determining 
the stability region of the dimer phase it is sufficient to find the minimum value of expression (|42(l at the edge of the 
domain and compare it to the energy of the dimer phase in the absence of charge ordering, Fd = —{1 + --0/4. 

We obtain that for WN < 21 the minimum of 142(1 at the edge of the domain is reached for Vi n = v ex = ±1/2. 
Since F(±l/2, ±1/2) = W/A — 1/2 > Fd for such parameters the dimer state is stable against the transition to the 
charge ordered state. 

For WN > 21, the minimum at the edge domain is reached at the points 

1 iV-2 1 
V °* =± Z N -l-£ > " m=± 2> (43) 

where the value of the energy is equal to 

I I W(WN -2I)(N -2) 



2N 2 2N(WN + 1 -W) 



(44) 



In particular, it follows from (|43jl that for N — > oo the minimum of energy at the edge of the domain corresponds just 
to the state mentioned above (v ex = Vm/2 = 1/4). 

The value of expression l|44|l can be larger or smaller than Fd- Therefore, there is a region of instability of the dimer 
phase in the space of parameters (c?, N). It is determined by the system of inequalities 

/ I- J W(WN -2I)(N -2) ^ () 
2N 4~ + 2N(WN + I-W) < 

WN > 21. (45) 

We analyze the system of inequalities numerically and find that for N < 10 the conditions (|45|l are not satisfied 
for any d. Consequently, in a multilayer system with even N < 10 the dimer state is stable against transition to 
the charge-ordered phase at any values of d. For N > 10 the inequalities l|45|l are satisfied in certain interval of d: 
d c i < d < d C 2- At such N and d there will be complete charge ordering of the inner layers and no phase coherence 
between those layers. The dependence of d& and d C 2 on N is shown in Fig. ^ Since limjv^oo d c i(N) = 0, for N — > oo 
the domain of values d for which the dimer state can exist is restricted by the condition d > 1A5£b- However, at large 
but finite N there exists a region of small d (d < £bV^/N) in which the formation of the dimer state is possible. 

To conclude this section we discuss shortly the role of tunneling. It was shown in Ref. |25| that the interlayer 
tunneling unfavours the charge ordering in multilayer systems. Based on a comparison of the energy of the system for 
v = and v = ±1/2 the authors of that paper arrived at the conclusion that in the case when the tunneling parameter 
t exceeds 0.05e 2 /e€s the charge ordering becomes energetically unfavorable for any values of the parameter d/£s ■ 
In the presence of tunneling between layers the phase coherence is no longer spontaneous, since tunneling leads to 
fixing of the phase of the order parameter. In that case the superfluid state with a constant phase gradient cannot 
be realized, but superfluidity of a "soliton" type can exis1ii£*2§i2£. But, according to the results of Ref. 0, already for 
t > O.Ole 2 /e£b the Landau criterium is not fulfilled at any superfluid velocities. Therefore, tunneling plays more of a 
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negative than a positive role for realization of the superfluid state in multilayer systems,. On the contrary, the use of 
multilayer systems with parameters lying in the stability region of the dimer state (Fig. ^) and with a small value of 
the tunneling is more promising. 



V. CONCLUSION 



We have shown that a state with spontaneous interlayer phase coherence can be realized in multilayer quantum Hall 
systems in a certain range of parameters. The phase coherence in such systems has a two-dimensional rather than a 
three-dimensional character; namely, the system separates into pairs of adjacent layers, and coherence is established 
only within each pair, and coherence is absent between electrons belonging to diffierent pairs. If the number of layers 
N < 10, then for an arbitrary distance between layers the state with spontaneous interlayer phase coherence is stable 
against transition to a charge-ordered phase. If the number of layers N > 10, the state with spontaneous phase 
coherence is stable only for sufficiently small or sufficiently large distances between layers. In the intermediate region 
of interlayer distances in a system with more than 10 layers a transition occurs to a charge-ordered state in which 
interlayer phase coherence is absent. 

Since the phase coherence has a two-dimensional character, a multilayer quantum Hall system can be considered 
as a solid-state analog of ultracold atomic Bose gases in optical superlattices (see, e.g., Ref. l28jh Among new effects 
that can expect in such systems we would mention the nondissipative interlayer drag of superfluid flow, in analogy 
with similar effect predicted for bilayer Bose systems22i22i2i. 
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